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Abstract. Integrable difference equations commonly have more low-order conservation laws than occur 
for nonintegrablc difference equations of similar complexity. We use this empirical observation to sift 
a large class of difference equations, in order to find candidates for integrability. It turns out that all 
such candidates have an equivalent affine form. These are tested by calculating their algebraic entropy. 
In this way, we have found several types of integrable equations, one of which seems to be entirely 
unrelated to any known discrete integrable system. We also list all single-tile conservation laws for the 
integrable equations in the above class. 



1. Introduction 

A quad-graph equation is a scalar difference equation for u(fc, Z), where (fc, I) £l? , which is of the form 

J^(A:,/,uoo,uio,woi,wii) = 0. (1) 

Here denotes uik -\- i,l + j) and we assume that J- depends on all four of these values. Various 
approaches have been used to discover quad-graph equations that arc integrable. Having developed 
the bilinear formalism for continuous integrable systems, Hirota discretized the bilinear operators for 
several known integrable systems, obtaining difference equations that had soliton solutions built-in [1-3]. 
By contrast, Capel et al. focused on discretizations of plane wave factors for the singular integral 
equations that are ubiquitous features of continuous integrable systems [4-6] . Whereas these approaches 
used discretizations of problems that were known to be integrable, Adler, Bobenko and Suris (ABS) 
dealt directly with quad-graph equations without reference to continuous systems. They obtained a 
classification of all integrable quad-graph equations that are consistent on a cube (and thus admit a Lax 
pair), subject to certain nondegeneracy conditions [7,8]. The idea that consistency on a cube is a sufficient 
condition for integrability was proposed independently by Nijhoff [9] and Bobenko and Suris [10]. 

To make further progress, we adopt a different strategy. There is a systematic method for construct- 
ing conservation laws of difference equations; this has been used to identify low-order conservation laws 
of many integrable quad-graph equations [11,12]. From this work, we observe that integrable differ- 
ence equations tend to have more low-order conservation laws than nonintegrable equations of similar 
complexity. Although this observation is purely empirical, we use it to sift a large class of quad-graph 
equations, in order to find equations that admit 'extra' conservation laws. (This approach is dual to 
that of Levi and Yamilov, who recently obtained some necessary conditions for the existence of higher 
symmetries - which again indicate integrability - for certain types of quad-graph equations [13]). Having 
obtained a shortlist of possible candidates for integrability, we test their algebraic entropy. 

Zero algebraic entropy is a signature of integrability [14-16]. This occurs for afhne linear quad-graph 
equations when an arbitrary set of initial conditions produces polynomial growth in degree as one moves 
away from the initial points (see §3 for details). Linear growth in degree implies that the quad-graph 
equation is linearizable; all known integrable quad-graph equations that are not linearizable exhibit 
quadratic growth. The calculation of algebraic entropy is a diagnostic test, rather than a constructive 
method. For instance, Hietarinta discovered a quad-graph equation that is consistent on a cube, but 
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does not appear in the ABS list [17]. A calculation of algebraic entropy showed that growth in degree 
for this quad-graph equation is linear; separately. Ramani et al. found a clever linearization [18]. 

In the next section, we determine conditions for the existence of extra conservation laws for a large 
class of quad-graph equations. Algebraic entropy is calculated in §3, and we find that most of the sifted 
quad-graph equations exhibit quadratic growth in degree. For completeness, we list the conservation 
laws in §4, before discussing our results and their consequences in §5. 

2. Classification of integrable cases via conservation laws 

In this section, we examine the conservation laws of equations of the form 

Ull = eiuoo + A{uio) ~ e2A{uoi). (2) 

Here each is either 1 or —1, and A is a nonlinear complex- valued function that is assumed to be 
'differentiable enough' (so that as many derivatives as needed are well-defined, at least locally). Equations 
in this class lack the D4 symmetry of the ABS equations. The class includes some known integrable 
equations (such as the Lattice KdV equation) and is simple enough for a complete classification of 
integrable cases to be possible. Henceforth, we use Aij to denote A{uij). Conservation laws on a single 
tile satisfy the determining equation 

F{k + 1, I, Mio, uj) - F{k, /, uoo, woi) + G(k, I + 1, uqi, t^) - G(fc, /, uoo, "lo) = 0, (3) 

where uj denotes the right-hand side of Q. We solve Q by deriving a sequence of its differential 
consequences, each of which eliminates at least one unknown function from the previous equation in 
the sequence (see [11] for a fuller explanation). This leads to an overdetermined system of functional- 
differential equations that can be solved completely. Specifically, we apply the commuting differential 
operators 

-Ci = <9io — eiA'iQdoo, £2 — <9oi + eie2^oi<9oo, 
(where dij denotes d/duij) to obtain 

(e2^'io^oi^oo+ei^io^oi)9oo^'(fc, h "00, uoi) + (e2^io^oi^oo-eie2^oi9io)9ooG(A:, I, uqo, uio) = 0. (4) 
Dividing by 62^10^01, then differentiating with respect to uqi yields the partial differential equation 




whose general solution is 

F{k, I, uqq, uqi) = /i(fc, I, eiuoa ~ 62^01) + f2{k, I, uqo) + faik, I, uqi). (6) 

Without loss of generality, set fs — O (absorbing the resulting trivial conservation law into /2 and G). 
Then ^ amounts to 

(doo - -^dw)dooG{k, I, uqo, mio) = -dQ^f2{k, I, uoo), 

whose general solution is 

G'(fc, /, Moo, tiio) = 9i{k, I, eiUQo + Aio) + g2{k, I, uio) - /2(fc, I, uqo). (7) 

At this stage, it is convenient to substitute ^ and ([7]) into the determining equation ([3]), using the 
difference equation ^ to eliminate Moo- This puts the determining equation in the form 

/i(^fc + 1, I, eiuio - £2^11) - fi{k, I, Ull - ^10) + f2{k + 1, I, uio) - f2{k, I + 1, uqi) 

+ gi(k, I + 1, eiuoi + An) - gilk, I, un + £2^01) + g2{k, I + 1, un) - 32 (fc, I, "lo) 0. (8) 
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Applying doidn to ([5]), wc obtain 

eiA[^g'l(k,l + l, eiUoi+Aii^ = ezAfji^'/ (/c, /, un + 63^01) , (9) 

where is the second derivative of g with respect to its third argument. This condition holds trivially 
if gi is linear in the third argument, which leads to two 'universal' conservation laws for which F and G 
are each linear in Aij . These arc 

Fi(k, I, uqo, uqi) = (Velei)''^' ^4(e2Woo - V^i^^oi), 

Gi{k, I, Uqo, uio) = (Veie2)''^'e2(e2Wio + ^00), (10) 

and 

F2{k, I, Uqo, Uqi) = ( - ^£162)''^' ^e2(e2Uoo + \/eie2^oi), 

G2(fc, ^, UOO, Wio) = ( - Vei<^2)''^'4(£2Mio + ^oo), (11) 

In order to find all functions A for which there are additional conservation laws on a tile, we now 
restrict attention to the case 5" 7^ oQ Dividing ([9]) by A'q^, then applying the operator doi — e2^oi9ii 
and rearranging the result, we obtain 

g'{'{k,l + l,eiUoi+An) _ A[,A'^, + e2{A'„,y A'{, 
g'{{k, Z + 1, eiUoi + Aii) yl[,iA'ii(ei -es^oi^'ii)' 

It is convenient to write A'^ = £?(^.y) = Bij, so that A^^- = BijB[^ (which is nonzero, as Aij is a 

nonlinear function of u^j) and = [BijY B'lj + Bij[B'^j)^. Then 

/ + l,emoi + An) ^ + £2^01^(1 
5"(A:, ; + 1, eiuoi + All) ei - £2-801-811 

Applying the operator 

doi 77— 9ll = Soi TTl £1 ■ 



(12) 



A[, "'dAoi 'dAn 

to (|12[) gives (after simplification) 

(l-eie2BoiSii)(B^,'i -eie2B;'i) -Boi(s;i)' + eie2Bii(B^i)' = 0. (13) 

This is the classifying equation that yields all functions A for which there exist conservation laws other 
than PH)) and (|lip . As stands, the functions Bqi and i?ii arc thoroughly entwined, but this can be 
resolved by one further differentiation, which yields the necessary condition 

(i3-/Sg' = 0. 

A simple calculation shows that (5^^)^ is a nonzero quadratic function of By ; substituting this into (jl3p 
and solving the resulting conditions gives 

{Bl^f = cl{Bl + 1) + (1 + eie2)c2B,j, ci, C2 G C. (14) 

(Here and henceforth, arbitrary constants arc denoted c or q.) This splits into four cases, as follows. 

Case I: ci = 0. 

In this case, we require £2 = ci and C2 7^ 0, in order that B[j is nonzero. Then 

K,^B,,^^{A,,+c:^, 



^If g" = but /" 7^ 0, similar calculations lead to precisely the same classifying equation II14I I. so nothing is lost by 
this assumption. 
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Aij = -. r-C3, where c = 2/c2 7^ 0. (15) 

(C4 - Uij) 

Then the sohition of ([T2|) . after absorbing the hnear terms mto /2 and .92, is 

5i(fc, Z + 1, eiiioi +^11) = C5hi(eiC4 - C3 - (eiiioi +^11)), C5 7^ 0. (16) 

This satisfies (O if ei = £2 = 1, but if ei = e2 = — 1 then ^ gives the further constraint C3 = 04. So this 
case leads to two possible equations, namely 

uii=uoo-c(— —), (17) 

\Uio Uqi/ 

and 

wii = -woo-c( \ ), (18) 

where the constant C4 has been absorbed into Uy . The point transformation uqq i-^ (— l)'^woo maps (|18p 
into (|17p. which is the lattice KdV equation (simplified slightly from the form stated in [19]). By solving 
([5]) for the remaining unknown functions, we obtain five conservation laws for (|17p . which are listed 
later (after all integrable quad-graph equations of the form © have been identified). The corresponding 
conservation laws for (|18p follow from the above transformation. 



Case II: ci 7^ 0, £2 = ei, = cf. 

In this case set C2 = £304, where £3 = ±1. Then the general solution of ([Till leads to the result 



1 - z"^ ' 



where the notation = gCiMij+c4 jg ugg^i henceforth. Then p2|l amounts to 

5i'(fc, Z + 1, eiuoi+^ii) £3-2^16-1^11+^3 
If £3 = £1, the general solution of (PO)) is 

a(fc, Z + 1) exp {ci(£iUoi + ^11) + C3 + £104} 



(19) 



(20) 



g'{(k, l + l, £iWoi+^ii) 



[1 - £iexp {ci(£iUoi + ^11) + C3 + £iC4}]' 



where a{k,l) is an arbitrary nonzero function. So when £3 = £1 = 1, the condition ^ gives only 
a{k, I) — a{k), whereas when £3 = £1 = —1 it also gives the constraint e^^'^3-c4) — i 'Writing ^ in terms 



of Zij, we obtain 



Zll _ Zqi - 1 

Zoo zio — 1 

and 



(when £3 = £1 = 1), (21) 



zoo^ii = 7 ^1 w"^ 7T , (when £3 = £1 = -1). (22) 

(2:10 — ^)(zoi — tj 

In [20] , Ramani et al. show that ([21]) is equivalent (under a point transformation) to the 'discrete Lotka- 
Volterra equation of type I' that was discovered by Hirota and Tsujimoto [21]. Levi and Yamilov recently 
found higher symmetries, a Lax pair and two conservation laws for a variant of this equation [13]. 
When £3 = — £1, equation (PU]) yields 

g'{{k, I + 1, £iUoi + An) = a{k,l + l)exp { - ci(£iUoi + An) - C3 - £104} 
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Equation Q produces the constraint a(fc, I) ~ a{k) when £3 = — ei 
gives a(fc, I) = 0. So wc obtain only one further equation, namely 

zii _ ziojzQi - 1) 
Zoo 2:01(210-1)' 
Note that (PT|) . (1^^ and are afhne linear in each 



-1, but when £3 
(when £3 ~ — ei = —1). 



1 it 



(23) 



Case III: ci 7^ 0, £2 = ei, 7^ cf. 

In this case 



Bij = C2 sinh [ciA.j + C3) - C2 



where C2 = C2/c^ and = 1 — 7^ 0. Then the general solution of A^j 



-By is 



and therefore ([T^ amounts to 
g'{'(k,l + l, eiuoi 



Cl 



C2(l-2.y) 

1 + eiC2 + eiC2exp {ci(eiUoi + An) + C3 + £104} 



g'{{k, l + l, eiUm+An) 



Hence 



and so 



g'{{k, l + l, eiuoi + All) = 



1 + eiC2 - eiC2Cxp {ci(eiMoi + An) + C3 + 6104} 
a{k, I + 1) exp {ci(eiUoi + An) + C3 + £104} 



1 + eiC2 - eiC2exp {ci(eiUoi + An) + 03 + eia] 
produces the constraint a{k,l) = a{k) when ei = 1. The resulting difference equation is 
£11 _ (210 + c)(2oi - 1) 
Zoo 



^{-1,0}, 



(24) 



(201 + C)(2l0 - 1) 

where c = (1 + C2)/(1 — C2). This is equivalent under a point transformation to the lattice MKdV 
equatioi0 (see [22-24]); in particular, a Lax pair for this equation is given in [24]. 

When ei = —1, the condition ([9|) yields a{k,l) = a{k), together with e'^i'^'i-'^a) — xhis leads to the 
difference equation 

(210 + c)(2oi + C) 



2002ll 



C I 



{-1,0}. 



(25) 



(210 - l)(2oi - 1) ' 

The point transformation 

zy ^ (-c)'^+'2(-i)''"" (26) 

maps (US]) to (gl). Note that when c = 0, (211) reduces to (gS]). Furthermore, is the hmit of ([M]) as 
c — *■ 00 with Zij fixed. So (|2ip and ((^5)) are each singular limits of the lattice MKdV equation. Moreover, 
the point transformation 

Zij ' ^ ^/zji 

maps (|23l) to the Lotka-Volterra type equation (PT|) FI 



(27) 



Case IV: ci 7^ 0, £2 = —^i- 

This is similar to Case III; the solution of (fTi)) is 

B^j = sinh (ciAy + C3). 

■^We thank Frank Nijhoff and Kenichi Maruno for alerting us to this. 
■^We are grateful to an anonymous referee for this observation. 
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Therefore 



and so 



Hence 



1 Zij 



g'('{k, I + 1, eiMpi + All) _ ci [l - ei exp {ci(eiMoi + An) + C3 + eiCj]] 
g'{{k, I + 1, eiuoi + An) [l + ei exp{ci(eiUoi + An) + C3 + £104}] 

a(fc, / + 1) exp {ci(eiUoi + An) + C3 + eic^] 



(28) 



g'l{k, l + l, eiUoi+^11) = 



1 + ei exp {ci(eiUoi + An) + C3 + 6104} 

— — 1 , and the resuhing difference equation 



(29) 



(30) 



When ei = 1. ([9]) gives the constraints a(fc, I) — a{k) and e 
is 

£11 _ _ (zio + l)(2:oi + 1) 
^00 (2^10 - l)(2oi - 1) ' 
When ei = —1, we obtain similarly a{k, I) = a{k), e^^* = —1, which leads to 

(zio + l)(zoi - 1) 

^00^11 — - 7 TTT — TT • 

(zio - f)(zoi + 1) 

Once again, the process has produced affine linear equations. It turns out that (|29p can be mapped 
to (|30|) by the point transformation with c = 1. 

3. Algebraic entropy 

To test the integrability of the previous lattice maps, we evaluate their algebraic entropy [25-27]. The 
system has an infinite dimensional space of initial conditions. We choose initial conditions on a diagonal 
regular staircase, which is shown in Figure I. 



A = {i 



e{0,I}}. 



(31) 



This defines a forward evolution towards the upper right corner of the lattice, and a backward evolution 
towards the lower left corner. 

The method is to let the system evolve, calculating Unm away from the diagonal by using (recursively) 
the defining relation on an elementary tile of the lattice. Each u„rn is a rational polynomial in terms 
of the initial conditions; the degree of the denominator is evaluated. The space of initial conditions is 
infinite-dimensional but, for any quad-graph equation, we need to specify only 2fc + 1 initial conditions 
to evaluate k iterates. This gives a sequence of degrees {dn}, as shown in Figure 1. The growth of that 
sequence gives the entropy 

e= lim - \n{d„). (32) 

n — ^00 n 

Vanishing of the entropy is the hallmark of integrability [14-16]. 

Although we are able to calculate only a limited number of terms of the sequence, it is possible to infer 
the exact value of the entropy. The reason is the existence of a finite recurrence with integer coefficients 
that is satisfied by the sequence of degrees. The most efficient way to find this recurrence is to fit the 
sequence with a Pade approximant. The existence of the recurrence on the degrees ensures that the 
generating function for the sequence of degrees is a rational fraction. 

Table 1 gives the sequences of degrees and the corresponding entropy for the various quad-graph 
equations in fJ21 Fot' comparison, we also include a nonintegrable equation that is only slightly different 
to P^ . namely 

(zio + 2)(zoi +2) 



200 Zll = 



2(zio-l)(2oi-l)' 



(33) 
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Figure 1. The distribution of degrees over the lattice. 



Equation 



Sequence 



{dn} 



m 

(ED 



(121 



1,3,7,13,21,31,43,57,.. 
1,2,4,7,11,16,22,29,.. 
1,3,6,10,14,18,22,26,.. 
1,3,6,11,18,27,38,51,.. 
1,3,7,13,21,31,43,57,.. 
1,3,7,13,21,31,43,57,.. 
1,3,7,17,41,99,239,577,. 



1 + n + n- 
1 + {n^ + n)/2 
4n- 2, (n > 2) 
+ 2, (n > 1) 
1 + n + n'^ 
1 + n + n? 
((l + \/2)"+i + (l-\/2)"+i)/2 



hi(l + ^/2) 



Table 1. The sequence of degrees for each equation. 
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Equation has linear growth of the degree, which indicates that this equation is hnearizable. This 
result is confirmed by the existence of an infinite family of conservation laws on a single tile (see §4). 
However, we have not been able to discover a linearizing transformation. The point transformation 
Zij 1-^ i/zij simplifies (|22|) to 

zoozii = [ziQ - l){zoi - 1), (34) 

but this is no more tractable than the original equation. 

All other cases that have more than two conservation laws exhibit quadratic growth of the degree, and 
therefore are claimed to be integrable, but not linearizable. This raises the question of whether any of 
the new integrable quad-graph equations can be mapped to any known equation. This will be discussed 
in §5. 



4. The conservation laws 

Although the lattice KdV equation (fT7|) , the lattice MKdV equation (p4|) and the discrete Lotka-Volterra 
equation (j2ip are not new, their conservation laws (on a single tile) have not previously been listed. 
Throughout this section, the universal conservation laws ([TU)) and pT|) span the first two conservation 
laws in each list (up to the addition of trivial conservation laws). So the 'extra' conservation laws are 
{Fi,Gi), where « > 3; here Fi denotes Fi(k, I, zqOi ^qi) and Gi denotes Giik, I, zqo, ^lo)- Lattice KdV 
has the following conservation laws. 

Fi = itoo + c/uoi, Gi = uio - c/uio; 

F2 = (-1)'+' (uoo + c/«oi) , G2 = (-1)^+'+! (uio + c/uw) ; 

F3 = In (uoo + c/uoi), G3 = In(wio); 

Fi = In(iioo), G4 = In (wio - c/uoo); 

F5 = -fcln(uoo) +/ln(uoo + c/uoi), G5 = -fcln(uio - c/uqo) + (l - l)ln(uio). 

For lattice MKdV, the single-tile conservation laws are: 

F, = In ( 'Jd^21^] , Gi = In ^ ''"^''^ + 

Fo 



zoi+c J' V ^10 - 1 

(-1)^+' In ( , G, =(-1)^+' In ( ^^oi^^o^l) 

\ZQa(zQi- I) J \ zio + c 



, I ZqqZoI - Zqq - Zoi - C\ 

F3 = In ■ , Ga = ln(zio + c); 

zoi+c ) 

, ,'2oo + c\ ^ , /^^oo^io + czoo + czio 

bA — In , G4 = In 



200 /' V zio(zoo + c) 



^ ( z,o{zo^+c) \ (zoozoi - zoo - 201 - c)2 



[zoi - l){zoo + c)^ J V 200(201 + c) (201 - 1) 
= fcln ( . zMzio-l){zoo+cr ^\ ^ / (2io + c)(2io-l) \ ^ f z,o 



(zio + c) (200210 + czoo + c2io - c)2y y J \{zio 
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The single-tile conservation laws for the discrete Lotka-Volterra equation (pij) are 
Fi = In (2:00(201 - 1)), 

F2 = (-l)'=+'ln(zoo(2oi - 1)), 
F3 = 200(1 - 201), 



= In(zoo), 



Gi 


= ln ( ■ 


G2 


= (-1) 


G3 


=210; 


Gi 


= ln ( ■ 


G5 


=fcln ( 



210 



200 + 210 -1 



F, = fcln ( ^) + nn (zoo(2o: - 1)), G, =fcln ( --10(210 - 1) \ ^ / £ic^\ ^ 

V201-I/ Vl^oo + 210 - 1)^/ V ^10 / 

Levi and Yamilov [13] recently derived an alternative form of ([21]) and listed two of its conservation laws, 
which are equivalent to (fi,Gi) and {F^.d). 

We now list the conservation laws corresponding to the remaining affine linear quad-graph equations 
that we have derived. Each of our equations that is not equivalent to Lattice KdV, Lattice MKdV or the 
Lotka-Volterra type equation is equivalent to either the linearizable equation (j22p or the new equation 
0. 

Equation (f22|) 

This equation has an infinite set of conservation laws, which depend upon two arbitrary functions a, (3: 



I? /; I iM / -'00201 - Zoo - 201 \ , /zoo2oi - 200 - z;oi \ ^ /, n 

Fa = a[l + 1) In + a(l) In , Gq = a[l) ln(l - zio); 

200 (,^01 ^ ^) / \ 201 



F, = m ln(l - zoi), G, = P{k + 1) In f f2£fHl_i2i_fH) + p^k) In ( »^i2£^ 

V 2oo(zio - 1) / V ^10 

This is a further indicator that, unlike the other quad-graph equations in our class, (|22p is linearizable. 
Equation dSO]) 

^1 - ( - 1)^'^^'"'^-^'-^^/^ In (^^) , Gi = cos {^) In (^-g^) + sin In (,^-^) ; 

. ( - m , G. = - sin (i-t^) m (^^) + cos In (,,^) ; 

^ 1^ (200 + 1)^(201 - 1) ^ ( (-1)' (200210 - 200 + 210 + 1)^(210 + 1) 



V 200(2:01 + 1) /' V 210(200 + 1)^(210-1) 



2oo(zoi + l)(2oi - 1) / ' V -210 



^ 200(201 + 1) \ ^j^^^^f (-l)''(2oo2oi + 200 - 201 + 1)^ 



(zoo + 1)^(201 - 1)/ V 200(201 + l)(2oi - 1) 

G. = fcln( , M):W2oo + m2io-l) ^ ^+nnr(-^" + ^)(-^"-^)Uh/ 



(210 + l)(2oo2lo - 200 + 210 + 1)2/ V 210 / V(^10 + l) 
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5. Comments 

Remarkably, although the original Ansatz ^ contained an arbitrary function A, each of the equations 
that we have found by sifting can be written in affine form, using a simple change of dependent variable. 
This has made the entropy calculation possible, because it gives rational evolution. 

It is natural to ask at this point how the equations that wc have derived compare to the known afhne 
linear quad-graph equations. Wc have already seen that in most cases, our Ansatz yields an equation that 
is equivalent under a point transformation to a known equation. Therefore it is important to characterize 
this equivalence, which can be done using the appoach introduced in [8]. Any affine linear quad-graph 
equation can be written in polynomial form: 

where w;, i = 1 . . . 4, are the values (of it.y or Zij as appropriate) at the four corners. For any choice of a 
pair of indices 1 < * < J < 4, define hij by 

hij{vk,vi) =dy,Q- dy^Q- Q ■ dy^dy.Q, j ^ k^l (35) 

It is then possible to associate to each of the four corners a polynomial 

rkivk) = {dv,h,,f - 2 h,, {dlh,,). (36) 

These polynomials play a central role in the classification of [8] , because (after a Mobius tranformation, 
if necessary), they can take one of six canonical forms, according to their root distribution. 
For example, the lattice MKdV equation yields (adjusting the notation for clarity) 

^^00^01 = (1 + c) (zio - 1) {zw + c) zii; 
^^00^10 = -(1 + c) (^01 - 1) (^01 + c) zii; 
^^00^11 = -(^01 - 1) (^;oi + c) (zio - 1) (zio + c); 

hzoizta = -(1 + cf Zoo Zii\ 

hzoizii ^ -(1 + c) (zio - 1) (zio + c) Zoo; 

hz^ozii = (1 + c) (^01 - 1) (^01 + c) Zoo. 

All of the functions hij are products of linear factors; this is the case for every equation in our classifi- 
cation. In other words, all of these equations are 'degenerate' in the sense used in [8]. Moreover 

roa = (1 + 0)" Zoo; 

rn = (l + c)4 zf,; 

no - {l + cf (l-zio)2 {zw + cf; 

roi = (l + cf (l-zoif {zoi+cf. 

These are in the canonical forms, but arc not in any of the cases that were classified in Theorem 2 
of [8]. Hence none of the equations that wc have studied are equivalent to any equation in the ABS 
classification. 

In summary, it is feasible to look for new intcgrable difference equations by searching for equations 
that admit 'extra' conservation laws. The class that wc have studied has been particularly fruitful, 
although only one of the equations (up to equivalence) seems to be unknown. A useful by-product is 
that one obtains a list of conservation laws, most of which are new (even for the known equations). The 
calculation of algebraic entropy is a clear indicator of integrability and linearizability. 

It is worth noting that our new equation (j30p is the only one with maximal asymmetry within the 
form of Ansatz ([2|), because ei = —€2 in this case alone. 

Two particularly important questions remain: does the new equation have a Lax pair description, 
and is it 3D-consistcnt? If wc wanted to check directly the consistency around the cube, we should first 
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choose an Ansatz for the form of the relations we want to use on the six faces of a cube. This leads 
one to ask which deformations of our models will be integrable. These might be Mobius transformations 
or other deformations which do not lie within the assumed Ansatz (0) ■ The analysis of the singularity 
pattern may be a way to tackle this problem. One should be prepared to accept deformed equations that 
are not affine; however, this is beyond the scope of the current paper. 
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